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Abstract 



A measurement of inclusive Z — > t + t production in pp collisions is presented, in 
the final states f/+hadrons, e+hadrons, e+}i, and p+fi. The data sample corresponds 
to an integrated luminosity of 36 pb _1 collected with the CMS detector at the LHC 
The measured cross section is a (pp — > ZX) x B (Z — >■ t+t - ) = 1.00 ± 0.05 (stat.) ± 
0.08 (syst.) ± 0.04 (lumi.)nb, which is in good agreement with the next-to-next-to- 
leading order QCD prediction and with previous measurements in the Z — > e + e~ and 
channels. The reconstruction efficiency for hadronic t decays is determined 
with a precision of 7%. 
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1 Introduction 



The measurement of the production cross section for pp — > ZX with Z — > t + t~ constitutes 
an important physics benchmark at the Large Hadron Collider (LHC). The T lepton can decay 
either into purely leptonic final states (t — > ev e v T denoted as "r e " or t — > }iv^v T denoted as 
"t^") or into hadronic final states denoted by "Th a d" consisting of a hadronic system and a v T . 
Constrained by the T mass, the hadronic system is characterized by a low particle multiplicity 
and a highly collimated jet which allows a Th a d signal to be separated from the large QCD jet 
backgrounds. The validation of the ih a d signal is essential in searches for new physics based 
on t leptons, such as Higgs boson decays to t + t~ [lj. The Compact Muon Solenoid (CMS) 
Collaboration recently reported a search for the Higgs boson in this channel |2J. Tau leptons 
can also be important signatures for searches of supersymmetry extra dimensions, and extra 
gauge bosons O. 

The Z — > t + t~ production cross section has been previously measured in proton-antiproton 
collisions by the CDF and DO Collaborations DUET m this study, Z — > t + t~ events in the T^Thad/ 
T e ihad/ fe^/ and t^t^ final states are selected from a sample of \/s = 7 TeV proton-proton 
collision data recorded by the CMS experiment at the LHC. The data sample corresponds to 
an integrated luminosity of 36 ± 1 pb _1 . The results are compared to previous measurements 
made in the e + e and final states [6|, providing a validation of T^d reconstruction and 

identification [7-9 1 and a direct measurement of the tau selection efficiency. 



The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal di- 
ameter, providing a magnetic field of 3.8 T. Within the field volume are the silicon pixel and 
strip tracker, the crystal electromagnetic calorimeter (ECAL) and the brass/ scintillator hadron 
calorimeter (HCAL). Muons are measured in gas-ionization detectors embedded in the steel 
return yoke. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the x axis pointing to the centre of the LHC, the y axis pointing up perpendicular to the LHC 
plane, and the z axis along the counterclockwise-beam direction. The polar angle 6 is measured 
from the positive z axis and the azimuthal angle (p is measured in the xy plane. Variables 
used in this article are the pseudorapidity t] = — ln[tan(0/2)], and the transverse momentum, 



The ECAL is designed to have both excellent energy resolution and high granularity, which 
are crucial for reconstructing electrons and photons produced in t decays. The ECAL is con- 
structed with projective lead tungstate crystals in two pseudorapidity regions: the barrel ( \rj \ < 
1.479) and the endcap (1.479 < \t]\ < 3). The transition regions between the barrel and the 
endcaps, 1.444 < \t]\ < 1.567, are not used for electron reconstruction. In the barrel region, 
the crystals are 25.8 Xo long, where Xo is the radiation length, and conform to a granularity 
of Arj x Acp = 0.0174 x 0.0174. The endcap region is instrumented with a lead-silicon-strip 
preshower detector consisting of two orthogonal strip detectors with a strip pitch of 1.9 mm. 
One plane is at a depth of 2Xo and the other at 3Xo. The ECAL has an energy resolution of 
better than 0.5% for unconverted photons with transverse energies above 100 GeV. The en- 
ergy resolution is 3% or better for the range of electron energies relevant for this analysis. The 
HCAL barrel and endcap regions cover the range 1 1] \ < 3 and are subdivided into towers with 
a segmentation of Arj x Acp = 0.087 x 0.087, corresponding to 5 x 5 ECAL crystals in the barrel 
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3 Lepton Reconstruction and Identification 



region. The HCAL forward region extends the calorimetry to \rj\ < 5. 

The inner tracker measures charged particle tracks within the range < 2.5. It consists of 
1 440 silicon pixel and 15 148 silicon strip detector modules and provides an impact parame- 
ter resolution of ~ 15 pm and a transverse momentum resolution of about 1.5% for 100 GeV 
particles. 

The muon barrel region is covered by drift tubes and the endcap regions by cathode strip cham- 
bers. In both regions resistive plate chambers provide additional coordinate and timing infor- 
mation. Muons can be reconstructed in the range < 2.4, with a typical pj resolution of 1% 
for p T w 40GeV. 

A more detailed description of CMS can be found in [ 10 1 . 



Muons produced by t decays in the Z — > t + t~ process are reconstructed in the tracker and 
muon chambers |[TT | . Quality cuts, based on the minimum number of hits in the silicon tracker, 
pixel detector, and muon chambers, are applied to suppress backgrounds from punch-throughs 
and decays in flight. 

Electrons are reconstructed by combining tracks produced by the Gaussian Sum Filter algo- 
rithm with ECAL clusters [T2| . Requirements are imposed that distinguish prompt electrons 
from charged pions mimicking electron signatures, and from electrons produced by photon 
conversions. 

The CMS particle flow (PF) algorithm [8| is used to form a mutually exclusive collection of 
reconstructed particles (muons, electrons, photons, and charged and neutral hadrons) by com- 
bining tracks and calorimeter clusters. These reconstructed particles are used to build compos- 
ite objects such as t's and jets, and to measure the missing transverse energy tfj. 

Electrons and muons from x decays are expected to be isolated in the detector, while leptons 
from heavy-flavour (c and b) decays and decays in flight are expected to be found inside 
jets. A measure of isolation is used to discriminate the signal from the QCD multijet back- 
ground, based on the charged hadrons, photons, and neutral hadrons falling within a cone 
AR = y 7 At] 2 + A(p 2 = 0.4 around the lepton momentum direction. In the T f/ T /( final state, a 
cone of AR = 0.3 is used. A sum of the pj for each particle type is made for charged hadrons 
with pj > 0.5 GeV and for photons and neutral hadrons with pj > 1 GeV; photons and neutral 
hadrons are excluded from the sum if they fall within inner cones of AR = 0.05 and AR = 0.08, 



respectively. The relative isolation variable is I^j 5 = E ( p T arge + p T on + py eutral ) / p T , where 



Vj i Pt ' an d Pj eu refer to the charged hadrons, photons, and neutral hadrons in the 
sum, respectively, and p T refers to the pj of the lepton £ = e, p. For muons, it is required that 
1*1 < 0.1, while for electrons it is required that 7 r p el < 0.08 in the barrel and 1^ < 0.04 in 
the endcaps. A similar formula is used for the T e r ;( final state but with the isolation quantities 
based directly on tracker and calorimeter information, calculated in a cone of AR = 0.3. In this 
case, it is required that 7 re i < 0.15 for muons and 7 re i < 0.1 for electrons. 

The Thad identification algorithm used in this measurement is known as the Hadrons Plus Strips 
algorithm [9], which starts from a high-px charged hadron and combines it with other nearby 
charged or neutral hadrons to reconstruct T decay modes. The identification of tz° mesons 
is enhanced by clustering electrons and photons in "strips" along the bending plane to take 
into account possible broadening of calorimeter signatures by photon conversions. To reduce 
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the contamination from QCD jets, the Thad-candidate isolation is calculated in a cone of AR = 
0.5 around the reconstructed r-momentum direction. It is requested that there be no charged 
hadrons with pj > 1.0 GeV and no photons with Ej > 1.5 GeV in the isolation cone, other than 
the t constituents. 

4 Event Selection 

The events preselected for this analysis in the T^Thad/ T e T }i, and t^t^ final states are required 
to pass the single-muon Level-1 (LI) trigger with a pj threshold of 7 GeV and the single muon 
High Level Trigger (HLT) [13], with a threshold varying from 9 GeV to 15 GeV, depending on 
the instantaneous luminosity. Events in the T e Th a d final state are selected using the single- 
electron LI trigger with a threshold of 8 GeV in the transverse energy, and a single-electron 
HLT trigger with a threshold of 12 GeV at higher instantaneous luminosity. At the end of the 
2010 data taking period, a combined e+r trigger was used in order to keep the rate low enough 
without a further increase of the electron threshold. The trigger has the same LI requirements 
as the electron trigger, but the HLT requires the presence of an electron of transverse energy 
larger than 12 GeV and a hadronic tau decay with pj larger than 15 GeV, tagged with a sim- 
plified version of the tau reconstruction algorithm with a less restrictive selection than the one 
used in the offline analysis. 

Offline event selection starts with the requirement of a well-defined primary vertex [14J. For 
the T^Thad and f" e Th a d ^ na ^ states, one isolated muon or electron is required with pj > 15 GeV 
and \rj \ < 2.1. The associated Th a d must be oppositely charged, with pj > 20 GeV and \rj\ < 2.3. 
In order to reject events coming from the W+jets background, the transverse mass of the lepton 

and the %i, Mj(£, St) = J 2p T #j ■ (1 — cos A(p) r is required to be less than 40 GeV, where A(p 
is the difference in azimuthal angle between the lepton and vectors. The Mj distribution 
and the selection requirement applied are illustrated in Fig. [I] (left) for the T f ,ih a d final state. 

For the T e T« final state, one isolated muon with \t]\ < 2.1 and one oppositely charged isolated 
electron with \tj\ < 2.4 are required, both with pj > 15 GeV. Further background suppression 
is achieved by requiring M T (}i, # T ) < 50 GeV and M T (e, # T ) < 50 GeV. 

For the r^Tp final state, events with two oppositely charged isolated muons with \rj\ < 2.1 are 
selected if one satisfies pj > 19 GeV and the other pj > 10 GeV. The requirement tfj < 50 GeV 
suppresses tt and W+jets backgrounds, and a requirement on the azimuthal angle difference 
between the muons, Acp^ > 2, rejects QCD background events, for muons originating from the 
same quarkonia decay or from a decay chain of heavy-flavour hadrons. Effective suppression 
of the Drell-Yan background is achieved with a multivariate likelihood ratio technique. For 
each event, the relative probabilities to belong to two event classes (Z — > t + t~ — > t^t^ and 
Z/7* — > are computed, exploiting the following set of discriminating variables: the 

ratio of the transverse momentum of the dimuon system to the scalar sum of the momenta 
of the two muons; the significance of the distance of closest approach (DCA) between the two 
muon helix tracks; the pseudorapidity of the dimuon system; and the azimuthal angle between 
the positive muon momentum and A normalized likelihood discriminant is computed and 
is shown in Fig. [l] (right). Events with a likelihood discriminant value larger than 0.87 are kept 
in the final sample. 

For the plots in Fig. [TJ the QCD multijet, W and diboson backgrounds were simulated with 
PYTHIA [15J, the Drell-Yan signal and background with the next-to-leading order (NLO) 
Monte Carlo generator POWHEG |T6til8| , and the top samples with Madgraph [19]. The tau 
decays were performed with Tauola pOl . The samples were normalized using the cross section 



4 



5 Background Estimation 
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Figure 1: Muon+^x transverse mass (left) for the T ?i Th a d final state. Likelihood discriminant 
(right) for the t^t ? , final state. The arrows show the final selection criteria applied. 



at next-to-next-to-leading order (NNLO) for Drell-Yan and W, at leading order (LO) for QCD, 
and NLO for the tt sample. The EWK label in the figure refers to the Z — > e + e~, }i + }i ~ , and 
diboson backgrounds. 

5 Background Estimation 

The Z — > t + t~ signal is established using the visible mass, which is the reconstructed mass 
of the ^Th a d system in the Tf T^ final states, and the dilepton invariant mass for the T e r f , and 
TftTp final states. Due to the presence of neutrinos in the final state, the Z — > t + t~ visible mass 
peak extends across the range 30-100 GeV, which is considerably broader than the Z resonance 
itself. The backgrounds generally span the same mass range, so an accurate determination of 
the signal yield requires effective background estimation techniques. The main background 
sources are QCD multijet processes, W+jets, and Z — > £ + £~, with small contributions from 
top-quark decays and dibosons. All backgrounds are measured in control regions where their 
contributions are enhanced and extrapolated to the signal selection region using selection effi- 
ciencies determined either from data or the Monte Carlo simulation. In the T^Thad final state, the 
QCD multijet background is estimated using samples of same-sign (SS) and opposite-sign (OS) 
events from data, for which the electron or muon isolation requirement is inverted; the QCD 
background estimate is based on the ratio of OS to SS yields. The visible mass distribution of 
the SS events is used to describe the background in the signal region. The W contribution is 
extracted from the region Mj(£, Ui) > 60 GeV, where it dominates the sample. In the t 6 t^ final 
state, the background contributions are expected to be small and are estimated from the simu- 
lation. The large-transverse-mass region is used to check the estimated background contribu- 
tions from W, diboson, and tt backgrounds. For the T„T« final state, the Drell-Yan muon-pair 
production events are selected with a reduced likelihood without including the muon DCA 
significance, and the resulting muon DCA significance distribution is fitted with signal and 
background shapes. The QCD multijet background estimate is obtained from a sample of SS 
dimuon events. 

The numbers of selected events and the expected background contributions are summarized in 
Table [TJ The statistical and systematic uncertainties in the methods used are also given, where 
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the uncertainties are added in quadrature. 



Table 1: Numbers of expected background events and number of data events passing all the 
selection criteria in the four final states. The uncertainties shown include the statistical and 
systematic uncertainties added in quadrature. 











T H T f> 

(M m < 70GeV) 


Z — > £ + £ , jet misidentified as T 


6.4 ± 2.4 


15.0 ± 6.2 




Z — > t + t~ , lepton misidentified as t 


12.9 ± 3.5 


109 ± 28 


2.4 ± 0.3 


20.1 ± 1.3 


tt 


6.0 ± 3.0 


2.6 ± 1.3 


7.1 ± 1.3 


0.15 ± 0.03 


W -> iv 


54.9 ± 4.8 


30.6 ± 3.1 


1.5 ± 0.5 


2.5 ± 2.5 


W — y tv 


14.7 ± 1.3 


7.0 ± 0.7 


QCD multijet 


132 ± 14 


181 ± 23 


ww/wz/zz 


1.6 ± 0.8 


0.8 ± 0.4 


3.0 ± 0.4 




Total background 


228 ± 16 


346 ± 37 


14.0 ± 1.8 


22.8 ± 2.8 


Total data 


517 


540 


101 


58 



6 Systematic Uncertainties 

The efficiencies for electron and muon reconstruction, identification, and isolation, as well as 
the trigger efficiencies are obtained from data. Correction factors for the values extracted from 
the simulation are determined using the method described in Ref. (tag-and-probe method). 
The measured efficiencies have a small dependence on pj and cover the full range of pj used 
in the analysis. The uncertainties on the correction factors are in the range of 0.2-1.1%. 

A similar technique is used to estimate the hadronic tau identification efficiency. A data sam- 
ple of taus is selected using Z — > t + t~ — > T /( Th a d events. The events are preselected without 
applying the full tau identification but only kinematic cuts and a set of requirements to sup- 
press the background from Z — > p~ + }i~ , W, and QCD events. The efficiency is then calculated 
as a ratio of the number of events that pass the tau identification requirement to the number 
of all preselected events. The total uncertainty of the measurement arises from the statistical 
uncertainty of the sample and the systematic uncertainties related to the understanding of the 
backgrounds in the preselection sample and amounts to 23% [9|. 

To estimate the efficiency of the Mj selection and the likelihood selection efficiency for the t^t ?( 
final state, an embedded sample is used where the muons in a Z — > data sample are 

replaced by simulated tau decays with the original muon momentum. The estimated uncer- 
tainties amount to 2%. 

To estimate the effect of the energy-scale uncertainties on the acceptance, the energy of all 
reconstructed objects (electrons, muons, taus, and jets) is varied within their respective un- 
certainty. After each independent shift, the missing transverse energy is recalculated and the 
event selection is repeated. The event yield is compared to the nominal value and the relative 
difference is quoted as the systematic uncertainty. The systematic uncertainties are in the range 
from 1% to 3.5%. 

To obtain the acceptance corrections, the D6T and Z2 PYTHIA tunes [21J were used in the 
simulation. The effect of the use of different tunes on the final extracted cross section is smaller 
than 1% and is not included in the systematic uncertainties. 
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7 Cross Section Measurement 



The main source of theoretical uncertainty in the calculation of the experimental acceptance 
comes from the parton distribution functions (PDFs). The central acceptance value is obtained 
with the CT10 PDF fl22l . The uncertainty is estimated using the error sets of the PDFs: CT10, 
MSTW2008NLO EH, and NNPDF2.0 and amounts to 2%. 

The experimental and theoretical uncertainties are summarized in Table [2] 

Table 2: Summary of the sources of systematic uncertainties and their estimated effect on the 
measured Z — > t + t~ cross section. 



Source 




T e Thad 






Trigger 


0.2% 


3% 


0.2% 


0.3% 


Lepton identification and isolation 


1.0% 


1.1% 


1% 


1% 


Thad identification 


23% 




Efficiency of Mj selection 


2% 




Likelihood selection efficiency 




2% 


Acceptance due to ih a d energy scale, 3% 


3.5% 




Acceptance due to e energy scale, 2% 




1.6% 


1.6% 




Acceptance due to ]i momentum scale, 1% 


1% 




1% 


2% 


Luminosity 


4% 


Parton distribution functions 


2% 



7 Cross Section Measurement 

The cross section is obtained for each final state with the following formula: 

N 

a(pp -> ZX) x B(Z -> t+t-) = , (1) 

where N is the number of extracted signal events, A is the acceptance of signal events, e is the 
signal selection efficiency, B is the branching fraction of the T-decay mode considered [25 1, and 
C is the integrated luminosity II26II . 

The visible mass distributions of the T^Thad/ TeThad/ T e"fy and T f ,T^ final states are shown in Fig.|2] 
To extract the signal, a fit is performed using the visible mass shapes from the simulation, ex- 
cept for the QCD multijet and Z — > backgrounds, which are obtained from data. For the 
simulation shapes, a variation of the electron and tau energy scales within their uncertainties is 
considered. The effect of the muon energy scale is negligible. The background normalizations 
correspond to the numbers listed in Table [l] and are allowed to vary within the estimated un- 
certainties. The background yields and signal shapes shown in Fig. [2] are those obtained from 
the fitting procedure. 

Table 3: Acceptance, selection efficiency and fraction of selected events outside the generator- 
level mass window for the four final states considered. 





T^Thad 


T e Thad 






Acceptance A 


0.13 


0.12 


0.074 


0.16 


Selection efficiency e 


0.37 


0.23 


0.55 


0.17 


Mass window correction / ou t 


0.03 


0.03 


0.02 


0.01 
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Figure 2: Visible mass distributions of the T ?/ Th a d (top left), T e Th a d (top right), r e T^ (bottom left), 
and T fl r fl (bottom right) final states. 



The acceptances were obtained with the NLO QCD program POWHEG in the Z — > t + t~ mass 
region 60 < M T + T - < 120 GeV. Table [3] shows the acceptances and the selection efficiencies 
for the different final states considered. The number of extracted events from the fit, Ng t , is 
corrected for the fraction of signal events outside the generator-level mass window, f out , where 
N = Not • (1 — /out) in Eq.[T] The correction factors used are also shown in Table [3] 

The measured values of the cross section from the four final states considered are shown in 
Table |4| where the uncertainties shown are due to statistical, systematic, integrated luminosity 
and t identification uncertainties. 

The measured values are compatible with each other and with the NNLO theoretical predic- 
tion, 0.972 ± 0.042 nb |27). They are also consistent with the CMS measurement based on 
Z —> e + e , events KH- 

The dominant uncertainty on the Z — > t + t~ cross section measurement comes from the Th a d 
reconstruction and identification efficiency. A simultaneous fit to all four final states is per- 
formed to obtain the cross section and a scale factor for the ih a d efficiency, which is the ratio 
of the efficiency in the data to that in the simulation. The result of the global fit is shown in 
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7 Cross Section Measurement 



Table 4: The measured values of the cross section from the four final states considered. The 
statistical, systematic and luminosity uncertainties are given. The uncertainty associated to the 
Thad reconstruction and identification efficiency, Th a d — ID, is shown separately. 



Final state 


a (pp -^ZX)xB(Z^ t+t") nb 


stat. 


syst. 


lumi. 


tID 




0.83 


0.07 


0.04 


0.03 


0.19 


felhad 


0.94 


0.11 


0.03 


0.04 


0.22 




0.99 


0.12 


0.06 


0.04 






1.14 


0.27 


0.04 


0.05 





Fig. [3] (left), where the likelihood contours for the best estimates of the cross section and the 
Thad efficiency scale factor are shown. In addition to the one-standard-deviation contour, the 
contours for which the likelihood L is reduced by 2A In L = 2.30 and 6.18 compared to the maxi- 
mum value are also shown, corresponding to a coverage of 68% and 95% in the two-parameter 
plane, respectively. The value of the cross section extracted from the fit is 

a (pp -^ZX)xB(Z^tV)= 1.00 ± 0.05 (stat.) ± 0.08 (syst.) ± 0.04 (lumi.) nb, 

which is compared to the individual final state measurements in Fig. [3] (right). The value for 
the cross section is dominated by the dilepton final states which have smaller systematic un- 
certainties. In the simultaneous fit, the Thad-ID scale factor is measured to be 0.93 ± 0.09. 

A more precise value of the hadronic tau reconstruction efficiency can be obtained by perform- 
ing a fit of the T f ,Th a d and T e Th a d final states, where the cross section is fixed to the value mea- 
sured by CMS in the electron and muon decay channels, a (pp — > ZX) x B (Z — > e + e , ]i +yr ) 
= 0.931 ± 0.026 (stat.) ± 0.023 (syst.) ± 0.102 (lumi.) nb 0. The extracted value of the T had -ID 
scale factor is 0.96 ± 0.07, which corresponds to a tau identification efficiency of (47.4 ± 3.3)% in 
data, for hadronically decaying tau leptons in the Z — > t + t~ sample with visible px > 20 GeV 
within the detector acceptance. 
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Figure 3: Likelihood contours for the joint parameter estimation of the cross section and the 
T-identification (left). The fitted central values (solid line) and their estimated lc uncertainties 
(dashed lines) are also shown. Summary of the measured Z — > t + t~ cross sections in the 
if/Thad/ TeThad/ T e T^, and T fi T f , final states, in the invariant mass range of 60 < M T + T - < 120 GeV 
(right). The inner error bar represents the statistical uncertainty. The extracted cross section 
from the combined fit and the NNLO theoretical prediction are also shown. 



9 



8 Conclusions 

A measurement of the cross section for the process pp — > ZX with Z — > t + t~ has been per- 
formed, based on the T f( Th a d/ T~ e Th a d, T e T^, and t^Tm final states. A clear signal is established 
in the visible mass distributions for all channels. The measured cross section, a (pp — > ZX) x 
B (Z -4 t+t - ) = 1.00 ± 0.05 (stat.) ± 0.08 (syst.) ± 0.04 (lumi.) nb, is consistent with theoret- 
ical expectations, and the CMS measurement for the Z — >■ e + e and decay channels. 
A global fit of the T^Thad and T e Th a d channels provides a 7% constraint on the efficiency for 
reconstructing hadronic tau decays in CMS. 
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